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Abstract 

A charge-symmetry-breaking nucleon-nucleon force due to the up-down 
quark mass difference is evaluated in the quark cluster model. It is 
applied to the shell-model calculation for the isovector mass shifts of 
isospin multiplets and the isospin-mixing matrix elements in IsOd-shell 
nuclei. We find that the contribution of the quark mass difference effect 
is large and agrees with experiment. This contribution may explain the 
Okamoto-Nolen-Schiffer anomaly, alternatively to the meson-mixing con- 
tribution, which is recently predicted to be reduced by the large off-shell 
correction. 

PACS numbers: 21.10.Hw, 24.85. +p, 21.60.Cs 
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Charge-symmetry-breaking (CSB) terms is required in the nuclear force for explaining 
several phenomena, e.g., the difference between the proton-proton and neutron-neutron 
scattering lengths [0 and the anomaly of mass differences of several mirror nuclei, called 
Okamoto-Nolen-Schiffer (ONS) anomaly [fj, [3], |4|, [|]. These experimental data imply that the 
nuclear force between two neutrons is slightly more attractive than between two protons. 

Three of the present authors have made an extensive analysis of isovector mass shifts and 
isospin-mixing matrix elements in IsOd-shell nuclei. It was shown that experimental values 
of these quantities are well explained by a short-range CSB force, but not by a long-range 
force. According to these findings, we look for the origin of a short-range CSB force in this 
report. 

The experimental isovector mass shift, b(u,T), is calculated by letting the isospin multi- 
plet mass equation, 

E(u, T, T z ) = a{u, T) + b(y, T)T Z + c(y, T)T 2 Z , 

reproduce the masses, E(u,T,T z ), of the 2T + 1 members of the nuclear isospin multiplet. 
In the above equation, v represents the quantum numbers other than the isospin, T, and its 
third component, T z . 

Several well-known CSB sources have contributions to the isovector mass shifts. Before 
considering the nuclear-CSB force, we subtract the contributions of the electromagnetic 
interactions (EM) and of the isovector single particle energy (ISPE). They are calculated by 
the analysis of IsOd-space shell model |6|, [/J. The EM contribution is evaluated in taking 
account of the Coulomb force between the protons with the charge form factor correction, 
the electromagnetic spin-orbit force and the magnetic spin-spin contact interaction with the 
magnetic form factor of the nucleon. The explicit formulae are given in ref. ||. The ISPE 
represents the CSB interactions between the valence nucleon and the 16 0-core. Subtracting 
them, we obtain the "experimental" data, shown in Fig. 1, which is to be compared with 
the nuclear-CSB contribution. 

We find a systematic behaviour, called "zigzag behaviour", i.e., that the "experimental" 
values of the isodoublets are reduced for A = An + 1 while for A = An + 3 they are enhanced. 
It is demonstrated in the previous studies || [F] that a short-range CSB force is the most 
probable source of the "zigzag behaviour" . 

A short-range nuclear-CSB force can be provided by the exchange of the mixed p-uo 
complex. It is extremely short-ranged, and explains at least half of the "zigzag behaviour" 
0. The meson-mixing potential is also used to explain other phenomena, e.g., the ONS 
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anomaly M, PL 110]. Goldman et al. argued, however, that an off-shell correction reduces 



the meson-mixing amplitude by a large factor. The correction is so large as to kill the meson- 



mixing contribution. Other calculations 12, [13|, [14[ also verify such an off-shell effect. It 



is, however, still controversial, because another analysis indicates the strong off-shell effect 
being inconsistent with the observed (^-dependence of p-7* coupling |Tj|. Further studies 
seem to be necessary. 

In the present article, we assume that the meson-mixing contribution is negligible, and 
examine another short-range CSB force due to the up-down quark mass difference, called 
quark effect (QE). Precisely speaking, it is a contribution of the quark mass difference in 
the gluon-exchange interaction between the valence quarks, while the meson-mixing may 
contain effects of the quark mass difference as well Such a quark CSB force has been 
used to explain the difference between the proton-proton and neutron-neutron scattering 
lengths and the ONS anomaly [fT7 |. In this study, we apply the quark CSB force to the 



shell model calculation for the mass shifts of nuclear isospin multiplets. 

The quark CSB potential is calculated in the nonrelativistic potential quark model of 
baryons fll6| . The model consists of the standard interaction hamiltonian, which contains a 
quark confining potential as well as one-gluon exchange interaction. The confining potential 
is assumed to be independent of isospin, i.e., flavour of the quarks fl], [16| . The one-gluon 
exchange interaction contains an isospin-symmetry-breaking term in the hyperfine contact 
interaction, 
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where Aj, <7j and are the colour SU(3) generator, Pauli spin matrix and mass of the 
constituent quark, q i} and ot s is the strong coupling constant. Other terms are estimated to 



have negligible contributions to the isospin-dependent nucleon-nucleon (NN) interaction [16]. 
This contact interaction yields a short-range NN force which has the range of the nucleon 
size. 

The two nucleon system is represented by a quark-cluster wave function composed by 



two three-quark clusters [16[]. The internal wave function of each nucleon is approximated 
by a Gaussian, and the internal variables are integrated out to obtain the NN potential from 
the quark-quark interaction. The obtained potential depends on the NN relative coordinate 
as well as the spin and isospin quantum numbers. 

It is found in the calculation of Chemtob and Yang (see Fig. 2 of ref.|fL6|) that the local 
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term of the hyperfme contact interaction is the leading term and represents the whole CSB 
interaction approximately. In this article, therefore, we deal only with the local term of the 
hyperfme contact interaction, and the other terms are neglected. 

In order to apply the potential to the shell-model calculation, we have to prepare the 
potentials for higher partial waves than the S-states. They are numerically calculated [ |l8|j . 
and we obtain a quark CSB potential, 



C*' = t/f ^ exp(-| i3 V; 2 ) I[r.(l) + r.(2)]^[l - ■ 9(2)1 (2) 



where is the nucleon size parameter of the three-quark cluster, and r 12 is the distance 
between the centers of the two clusters, t x (i) is the third component of the Pauli matrix for 
the nucleon isospin, and a(i) is for the nucleon spin. 5m = m^ — m n is the up-down quark 
mass difference and m = \{m& + m u ) is the average of the masses. The input parameters 
are taken from ref. [16|, 5m = 6 MeV, m = 330 MeV, a s = 1.624 and j3~ l = 0.616 fm. 
The nuclear matrix element of the QE potential gives the isovector mass shift, 

b QE (^T) = 1 (v,T\\v QE \\v,T), 

<(2T + 1)T(T+ 1) 



in the first order perturbation, where the matrix element (v, T\\vqe\\u, T) is reduced with 
respect to the isospin. 

The nuclear wave function, \v,T), is calculated || with Wildenthal's effective hamil- 
tonian [19| in the complete IsOd-shell space. Because the short-range QE potential is inte- 
grated, the calculation is sensitive to the short-range structure of the relative wave function 
of the two nucleons. We include a short-range correlation by a correlation function of ref. pt)| 
multiplied to the relative two-nucleon wave function. 



The calculated QE contributions are shown in Fig. 2. The QE contributions are around 
100 keV, and the average ratio of QE/Coulomb contributions for 143 multiplets is 5.3%. This 
result is consistent with the expected contribution which is introduced phenomenologically 
to explain the ONS anomaly in literatures [Q, |5| , and consistent with the previous calculation 
of the quark effect for the anomaly ||17|| . 

The "zigzag behaviour" seen in the "experimental" mass shifts (Fig. 1) is also reproduced 
in the calculation (Fig. 2). Namely, the quark contributions are larger for the A = An + 3 
isodoublets, and smaller for A = An + 1. This behaviour is due to the short-range nature of 
the quark CSB force 0. 
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We need to comment on the calculation of the ISPE contribution. The ISPE used in Fig. 
1 is determined by a x 2 -fitting to the experimental mass shifts, including the QE contribution 
by the prescription of ref.@. By this procedure, the ISPE could be adjusted to produce the 
"zigzag behaviour" given by the QE contribution, even if the original data would not have a 
"zigzag" feature. To examine this, we have made another x 2 -fitting without any nuclear-CSB 
contribution, and still found the "experimental zigzag behaviour" (see ref.|| 0). 



The quark CSB potential of eq.(H) is also applied to off-diagonal matrix elements, i.e., 
isosp in-mixing matrix elements. The results are given in Table 1, with experimental values 
extracted from the strength of the isospin-forbidden beta decays. We find that the calculated 
matrix elements agree with most of the experimental values. 

To the off-diagonal matrix elements, the quark contribution is found to be large, compa- 
rable with the electromagnetic contributions (EM) in the table, which are roughly equal to 
the contributions of the Coulomb force. By contrast, the mass shifts are dominated by the 
Coulomb contributions. Such large contributions to the off-diagonal matrix elements can be 



attributed to the short-range feature of the quark CSB force plj, g. 

In particular, the beta decay of 24 Al (Table 1) indicates a large isospin-mixing in an 
excited state of the daughter nucleus 24 Mg. Our calculation also gives a large QE contribution 
to this matrix element, and thus is consistent with experiment. 



In summary, the isovector mass shifts and isospin-mixing matrix elements relevant to 
the isospin-forbidden beta decays are calculated considering the effect of the up-down quark 
mass difference in the direct gluon-exchange process. The charge-symmetry-breaking (CSB) 
nucleon-nucleon potential due to the quark mass difference is constructed with the quark 
cluster model, and is applied to the shell model calculation of the IsOd-shell space. 

The quark contributions are found to be about 5% of the Coulomb contributions to the 
isovector mass shifts of the nuclear isospin multiplets. It is consistent with the expected 
contributions & H to explain the anomaly of mass differences of the mirror nuclei, known 
as the Okamoto-Nolen-Schiffer anomaly. Also, the quark contribution is found to have a 
systematic behaviour in the contributions to the mass shifts of the isospin doublets, i.e., for 
A = An + 1 the mass shifts are reduced while for A = An + 3 they are enhanced. This 
behaviour is consistent with the experimental one. 

The calculated values of the isospin-mixing matrix elements are consistent with the most 
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of the experimental values extracted from the isospin-forbidden beta decays. The quark CSB 
force is found to have large contributions to the isospin-mixing matrix elements, comparable 
with the Coulomb contributions. A particularly large value of the experimental matrix 
element of 24 Mg is found to be explained by the large quark contribution. These findings 
strongly suggest the existence of a short-range CSB interaction. 

The quark mass difference effect is, of course, not a unique possible source of the short- 
range CSB force. However, it seems to have large contribution to the observables, and have 
a favourable feature for explaining the known data. 

The numerical calculations were performed with the FACOM M780 computer system 
at the Institute for Nuclear Study, University of Tokyo. This study is partly supported 
by a Grant-in-Aid for Scientific Research (05243204 and 06234206) from the Ministry of 
Education, Science and Culture (Monbusho), and one of the authors (S.N.) is supported by 
a scholarship from the Soryushi Shogakukai. 
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Table 1 Isospin-mixing matrix elements in keV and their decompositions into the contribu- 
tions of quark effect (QE), electromagnetic interactions (EM) and others. The EM consists 
of a large Coulomb contribution plus the other small contributions, and contains the isoten- 
sor terms of the electromagnetic interactions. "Other" consists of the ISPE contribution 



and the isotensor contribution due to the pion mass difference |25|| . Experimental values 
are extracted from the strengths of isospin-forbidden beta decays of the given initial states 

M-M- 



initial state 


19Q 


20 p 


24 A1 m 


24 Na 


24 Al 


27 Mg 


28 Mg 


QE 


-12.2 


11.5 


7.4 


5.5 


36.8 


2.6 


-8.2 


EM 


3.3 


8.3 


10.0 


6.2 


43.3 


-16.0 


-19.2 


other 


25.5 


-11.6 


-8.3 


-3.3 


-21.6 


41.1 


67.2 


total 


16.6 


8.2 


9.1 


8.5 


58.5 


27.8 


39.7 


experiment 


20(10) 


14 +29 


49(5) 


5.4(22) 


106(40) 


<J-°-3.6 


20.6(16) 



Figure captions: 



Fig. 1 The experimental isovector mass shifts p2| , [23] , |24 
after subtracting the electromagnetic contributions and 
the isovector single-particle-energy contribution in keV. 
The isodoublets (T = 1/2) are denoted by the filled cir- 
cles, and the others (T > 1/2) are by the crosses. 

Fig. 2 The quark contribution to the isovector mass shifts 
in keV. Notations are the same as in Fig. 1. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9502013vl 
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